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ABSTRACT 


An  improved  World  Geodetic  System  has  been  developed  to  satisfy 
the  mapping,  charting,  and  geodetic  needs  of  the  Department  of 
Defense.  The  system,  designated  WGS  72,  represents  the  culmination 
of  approximately  five  years  of  data  collection  within  the  Department 
of  Defense  and  the  scientific  community.  Principal  agencies  involved 
in  the  developmental  cycle  include  the  USAF,  Defense  Mapping  Agency, 

US  Naval  Weapons  Laboratory,  and  Naval  Oceanographic  Office.  The 
method  of  solution  for  an  adjustment  of  this  magnitude  is  characterized 
by  the  formation  of  a  large  scale  matrix  by  combining  normal  matrices 
from  each  of  the  major  data  input  sets.  The  scaling  and  weighting 
processes  for  the  final  matrix  are  discussed.  The  resultant  ellipsoid 
parameters,  datum  shifts  and  their  related  accuracies  are  presented 
and  compared  against  similar  quantities  in  other  recent  geocentric 
systems. 


THE  DEPARTMENT  OF  DEFENSE 
WORLD  GEODETIC  SYSTEM  1972 


INTRODUCTION 

The  Department  of  Defense  (DOD)  in  the  late  1950s  generated  a 
geocentric  reference  system  to  which  different  geodetic  networks  could  be 
referred  and  compatibility  established  between  the  coordinates  of  sites 
of  interest.  Efforts  of  the  Army,  Navy  and  Air  Force  were  combined  leading 
to  the  development  of  the  DOD  World  Geodetic  System  1960  (WGS  60).  In 
accomplishing  WGS  60,  a  combination  of  available  surface  gravity  data, 
astrogeodetic  data  and  results  from  HIRAN  and  Canadian  SHORAN  surveys  were 
used  to  obtain  a  best-fitting  ellipsoid  for  the  major  datum  areas.  The 
sole  contribution  of  satellite  data  to  the  development  of  WGS  60  was  the 
value  for  the  ellipsoid  flattening  (1/298.3  ±0.1),  which  was  obtained  from 
the  nodal  motion  of  Satellite  1958  6.  The  semimajor  axis  of  the  WGS  60 
Ellipsoid  was  determined  as  6,378,165  +50  meters. 

In  January  1966,  a  World  Geodetic  System  Committee  was  charged  with 
the  responsibility  of  developing  an  improved  WGS  needed  to  satisfy  mapping, 
charting  and  geodetic  requirements.  Additional  surface  gravity  observations, 
results  from  the  extension  of  triangulation  and  trilateration  networks  and 
large  amounts  of  Doppler  and  optical  satellite  data  had  become  available 
since  the  development  of  WGS  60.  Using  the  additional  data  and  improved 
techniques,  the  WGS  Committee  produced  WGS  66  which  has  served  DOD  needs 
since  its  implementation  in  1967.  The  defining  parameters  of  the  WGS  66 
Ellipsoid  were  the  flattening  (1/298.25  +0.02),  determined  from  satellite 
data  and  the  semimajor  axis  (6,378,145  +20m),  determined  from  a  combination 
of  Doppler  satellite  and  astrogeodetic  data  involving  a  geoid-match 
technique.  A  24th  degree  and  order  geopotential  coefficient  set  derived 
from  a  harmonic  analysis  of  a  worldwide  5°x5°  mean  free  air  gravity  anomaly 
field  was  selected  as  the  WGS  66  Gravitational  Model.  This  geopotential 
coefficient  set  was  also  used  in  a  spherical  harmonic  expansion  to  obtain 
the  Worldwide  WGS  66  Geoid.  Also,  a  geoid  referenced  to  the  WGS  66 
Ellipsoid,  providing  a  detailed  representation  for  limited  land  areas, 
was  derived  from  available  astrogeodetic  data.  Datum  shift  constants 
for  the  North  American  Datum  1927  (NAD  27),  European  Datum  (ED)  and 
Tokyo  Datum  (TD)  were  obtained  for  each  datum. 

The  same  WGS  Committee  began  work  in  1970  to  develop  a  replacement 
for  WGS  66.  Since  the  development  of  WGS  66,  large  quantities  of 
additional  data  had  become  available  from  both  Doppler  and  optical 
satellites,  surface  gravity  surveys,  triangulation  and  trilateration 
surveys,  high  precision  traverses  and  astronomic  surveys.  In  addition, 
greater  capabilities  had  been  developed  in  both  computers  and  computer 
software.  Further,  continued  research  in  improved  computational  procedures 
and  error  analyses  had  produced  better  methods  and  a  greater  facility  for 
handling  and  combining  data. 
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After  an  extensive  effort  extending  over  a  period  of  approximately 
three  years,  the  Committee  completed  the  development  of  the  Department  of 
Defense  World  Geodetic  System  1972  (WGS  72).  Selected  satellite,  surface 
gravity  and  astroqeodetic  data  available  through  1972  from  both  DOD  and 
non-DOD  sources  was  used  in  a  Unified  WGS  Solution  (a  large  scale  least 
squares  adjustment).  The  results  of  the  adjustment  consist  of  corrections 
to  initial  station  coordinates  and  qeopotential  coefficients.  This  paper 
presents  a  brief  review  of  the  data  and  methods  used  and  a  summary  of 
the  results. 


DATA  USED  IN  WGS  72  DEVELOPMENT 

The  largest  collection  of  data  ever  used  for  WGS  purposes  was 
assembled,  processed  and  applied  in  the  development  of  WGS  72.  Both 
optical  and  electronic  satellite  data  was  used.  The  electronic  satellite 
data  consisted,  in  part,  of  Doppler  data  provided  by  US  Navy  and  cooperating 
non-DOD  satellite  tracking  stations  established  in  support  of  the  Navy's 
Navigational  Satellite  System  [1].  Doppler  data  was  also  available  from 
the  numerous  sites  established  by  Geoceivers  (geodetic  receiver)  during 
the  1971-1972  time  period  [2].  Additional  electronic  satellite  data  was 
provided  by  the  SECOR  (Sequential  Collation  of  Range)  Equatorial  Network 
completed  by  the  US  Army  in  1970  [TJ.  Optical  satellite  data  from  the 
Worldwide  Geometric  Satellite  Triangulation  Program  was  provided  by  the 
BC-4  camera  system  [4]  [5],  The  Paker-Nunn  camera  data  of  the  Smithsonian 
Astrophysical  Observatory  (SAO)  was  also  used  [6]. 

Doppl er 


Geopotential  coefficients  and  some  station  coordinates  have  been 
determined  simultaneously  in  general  geodetic  solutions  using  sampled 
Doppler  observations,  taken  between  1962  and  1970,  of  satellites  at 
various  incl inations.  However,  an  accurate  (and  most  complete)  set  of 
station  coordinates  has  been  derived  by  the  point  positioning  method 
using  sampled  or  integrated  Doppler  observations  of  the  Navy  Navigation 
Satellites  (NAVSATs),  whose  polar  orbits  are  described  by  current 
mathematical  models. 

Since  1964,  locations  of  sites  occupied  by  mobile  Doppler  vans  have  been 
determined  using  a  minimum  of  40  passes  of  polar  NAVSAT  data.  Although 
coordinates  of  many  of  the  sites  have  been  subsequently  determined  together 
in  a  general  geodetic  solution  usinq  observations  of  several  satellites, 
an  independent  solution,  if  determined  by  current  point  positioning 
methods,  is  considered  more  accurate  because  of  the  sufficiency  of  data, 
the  high  accuracy  of  polar  orbit  computations,  and  the  precision  of 
NAVSAT  time  calibration.  With  this  in  mind,  the  semi- permanent  TRANET 
stations  were  independently  repositioned  in  1970  usinq  four  sufficient 
se'i  of  modern  data  in  a  way  that  was  consistent  with  the  Bureau 
International  de  l'Houro  (Bill)  pole  computations.  The  resulting  station 


2 


- 


coordinates,  defined  as  the  NWL-9D  geocentric  system,  is  the  basic 
Doppler  system  contributing  to  WGS  72. 

Normal  equations,  representing  satellite  data  from  approximately 
30,000  Doppler  and  500  optical  passes,  were  provided  as  a  matrix  and 
vector  of  dimension  822  allowing  corrections  to  unnormalized  gravitational 
field  coefficients.  Love's  constant,  and  the  geocentric  polar  coordinates 
of  114  observing  stations.  The  gravitational  field  coefficients  include 
all  coefficients  through  degree  and  order  19  (except  the  first  degree), 
zonal s  through  degree  24  and  additional  resonance  terms  through  order  27. 
Correlated  variations  in  orbit  and  bias  parameters,  implicit  in  the  com¬ 
putations  but  not  appearing  in  the  matrix,  are  fully  accounted  for  in  any 
solution. 

Normal  equations  for  stations  and  the  full  set  of  479  gravitational 
field  coefficients  were  developed  for  about  70  percent  of  the  satellite 
data  processed,  including  all  330  days  of  non-polar  and  92  days  of  polar 
Doppler  observations  and  all  96  days  of  optical  observations.  In  an 
effort  to  include  a  larger  number  of  Doppler  National  Geodetic  Satellite 
Program  (NGSP)  sites  in  the  matrix,  the  normal  equations  for  the  remaining 
30  percent  of  the  data  consisting  of  82  two-day  polar  arcs,  were  developed 
only  for  station  coordinates  and  three  gravitational  field  coefficients: 

Cq  q;  q;  C3  q.  Use  of  the  NWL-9B  gravitational  coefficients  set  insur¬ 
ed  an  RMS  orbit  error  of  only  three  meters  for  each  two-day  polar  arc, 
allowing  the  neglect  of  additional  gravity  variations  on  station  deter¬ 
mination. 

Stations  which  contributed  data  fall  into  the  semi -permanent  and 
mobile  classifications  (Figure  1).  Nearly  85  percent  of  the  data  was 
received  at  18  stations  of  the  semi -permanent  Doppler  network.  The 
remaining  15  percent  was  taken  by  six  similarly  equipped  mobile  vans 
operating  for  six  to  twelve  weeks  at  over  120  worldwide  locations.  Some 
20  of  these  sites  do  not  appear  in  the  matrix  because  of  size  limitations 
in  the  processing  program.  Additionally,  a  small  number  of  observations 
were  contributed  by  the  twelve  original  Baker-Nunn  stations  operated  by 
SA0. 


BC-4 


In  the  worldwide  network  observed  with  BC-4  cameras,  the  measurements 
of  approximately  2350  plates  were  used  to  establish  the  relative  positions 
of  49  tracking  stations  (Figure  2).  The  BC-4  camera  observations  were 
introduced  in  the  Unified  WGS  72  Solution  in  the  form  of  normal  equations 
developed  by  National  Ocean  Survey  (NOS)  for  their  TRI  Solution.  The 
procedure  used  by  NOS  was  to  develop  normal  equations  separately  for  each 
event;  i.e.,  a  set  of  two  or  more  simultaneous  photographs  of  the  satellite. 
The  unknowns  in  the  normal  equations  were  corrections  to  the  rectangular 
coordinates  of  the  satellite  and  the  observing  stations.  Normally,  there 
were  seven  satellite  positions  per  event  and  as  many  as  four  observing 
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Figure  1.  Doppler  Satellite  Ground  Stations  Providing  Data  for  WGS  72  Development. 


Worldwide  Geometric  Satellite  Triangulation  Network,  BC-4  Cameras 


stations.  The  satellite  position  unknowns  were  then  eliminated  leaving 
a  set  of  event  normal  equations  relating  the  coordinates  of  the  stations 
that  observed  that  event.  The  event  normal  equations  were  then  summed 
over  all  events  and  the  resulting  set  of  equations  used  in  the  Unified 
WGS  72  Solution.  A  total  of  856  two-station,  194  three-station  and  14 
four-station  events  were  used  to  relate  the  49  observing  stations. 

Baker-Nunn 


Several  geodetic  solutions  have  been  accomplished  by  the  Smithsonian 
Astrophysical  Observatory  using  data  obtained  from  the  optical  observations 
of  satellites  by  Baker-Nunn  cameras.  A  recent  solution,  the  1969  Smithsonian 
Standard  Earth  II  [7],  utilized  optical  data  from  Baker-Nunn  and  Mini¬ 
track  Optical  Tracking  System  (MOTS)  cameras  (Figure  3).  This  SAO  solution 
for  117  station  coordinates  (39  stations)  and  296  tesseral  harmonic  coefficients 
is  based  on  optical  and  laser  satellite  data  as  well  as  surface  mean  free 
air  g.avity  anomalies.  In  addition,  some  observations  of  deep-space 
probes  are  included.  The  complete  solution  represents  a  combination  of 
dynamic  and  geometric  satellite  data  with  surface  gravity  observations. 

The  normal  equations  developed  by  SAO  while  deriving  Standard  Earth  II 
were  incorporated  into  the  Unified  WGS  72  Solution.  Standard  Earth  II 
is  itself  the  result  of  a  combined  solution  that  incorporated  laser  and 
optical  satellite  tracking  data,  surface  gravity  observations  and  space 
probe  data. 

SECOR 


The  SECOR  measurements  used  ’n  the  Unified  WGS  Solution  for  WGS  72 
were  those  obtained  at  stations  forming  the  SECOR  Equatorial  Network 
(Figure  4).  The  results  of  the  SECOR  Equatorial  Network  Project  were 
incorporated  into  the  Unified  WGS  72  Solution  in  the  form  of  normal 
equations  resulting  from  short  arc  adjustments. 

Surface  Gravity 

The  surface  gravity  field  used  in  the  Unified  WGS  Solution  consisted 
of  a  set  of  410  10°xl0°  equal  area  mean  free  air  gravity  anomalies  deter¬ 
mined  solely  from  terrestrial  data.  This  gravity  field  includes  mean 
anomaly  values  compiled  directly  from  observed  gravity  data  wherever  the 
latter  was  available  in  sufficient  quantity.  The  values  for  areas  of 
sparse  or  no  observational  data  were  developed  from  geophysically  com¬ 
patible  gravity  approximations  using  gravity-geophysical  correlation 
techniques.  Approximately  45  percent  of  the  410  mean  free  air  gravity 
anomaly  values  were  determined  directly  from  observed  gravity  data.  Using 
the  surface  gravity  data,  normal  equations  were  formed  for  a  total  of  319 
gravitational  parameters.  These  included  the  zonal  harmonic  coefficients 
through  degree  20  and  the  tesseral  harmonic  coefficients  through  degree 
17  and  order  15.  Table  1  shows  the  quantity  of  gravity  data  available 
through  November  1972  for  WGS  development  and  evaluation  purposes. 
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i.  Baker- Nunn  and  MOTS  Stations. 


Figure  4.  SECOR  Equatorial  Network, 


Table  1 


Mean  Free  Air  Gravity  Anomaly  Information  -  November  1972  Field 


1' 

3x1°  Ag  Data 

5 

°x5°  &g  Data 

Region 

Number 

of 

Values 

RMS 

Value* 

(Mqals) 

Average  1 
Val ue* 

( Mga 1 s ) 

Number 

of 

Values 

RMS 

Value* 

(Mgals) 

Average 

Value* 

(Mgals) 

North  America 

10,557 

±26.9 

-1.4 

489 

±17.9 

-1.7 

South  America 

2,830 

36.0 

-1.4 

144 

18.9 

-0.5 

Europe 

5.256 

32.6 

0.0 

251 

20.2 

-1.4 

Asia 

7,871 

34.5 

-0.1 

353 

18.4 

-0.7 

Eurasia 

13,127 

33.9 

-0.1 

604 

19.1 

-1.0 

Australia 

1,618 

27.4 

63 

20.4 

-4.8 

Africa 

963 

21.1 

193 

15.4 

4.8 

Antarctica 

296 

39.9 

-4.7 

49 

22.1 

-8.7 

*AU  Ag  data  referenced  to  the  WGS  72  Ellipsoidal  Gravity  Formula. 
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Astrogeodetic  Geo i d  ! iei yht  Data 

The  astrogeodetic  data  in  its  basic  form  consists  of  deflection  of 
the  vertical  components  referred  to  the  various  national  geodetic  datums. 
These  deflection  values  were  integrated  into  astrogeodetic  geoid  charts 
referred  to  these  national  datums  using  the  methods  described  in  [8]. 

The  increase  in  the  holdings  of  astrogeodetic  data  since  1966  is  summarized 
in  three  reports  to  the  International  Union  of  Geodesy  and  Geophysics 
(IUGG)  Assemblies  in  1967  [8]  and  1971  [9]  [10].  The  astrogeodetic  geoid 
height  data  was  inserted  into  the  Unified  WGS  Solution  in  the  form  of  a 
430x430  normal  equation  matrix,  the  right  hand  side  being  a  430x1  column 
of  constants.  The  astrogeodetic  geoid  height  matrix  contributed  to  the 
Unified  WGS  Solution  by  virtue  of  providing  additional  and  more  detailed 
data  for  land  areas. 

Survey 

Normal  equations  representing  conventional  ground  survey  data  were 
included  in  the  Unified  WGS  Solution  to  enforce  a  consistent  adjustment 
of  the  coordinates  of  neighboring  observation  sites  of  the  BC-4,  SECOR, 
Doppler,  and  Baker-Nunn  systems.  Also,  normal  equations  representing 
eight  geodimeter  long  line  precise  traverses  were  included  for  the  purpose 
of  controlling  the  scale  of  the  unified  solution. 

To  relate  each  set  of  two  to  five  stations  located  up  to  100  kilometers 
apart,  a  normal  equation  submatrix  was  produced,  treating  the  survey 
coordinates  as  observables  with  a  given  standard  error  (usually  0.5  meter) 
and  allowing  as  parameters  the  coordinates  of  all  involved  stations  and 
seven  transformation  constants  for  the  survey  data:  three  rotation,  three 
translation,  and  one  scale.  Rotation  and  scale  parameters  were  given  non¬ 
zero  a  priori  standard  errors  (usually  50  ppm)  to  weaken  the  tie  over 
long  distances.  The  final  tie  matrix  used  was  the  sum  of  54  such  sub¬ 
matrices  after  individual  elimination  of  the  transformation  parameters. 

A  matrix  representing  the  eight  geodimeter  scalars  was  formed  for  coordi¬ 
nates  of  paired  BC-4  stations  located  at  endpoints  of  the  long  survey 
lines. 


UNIFIED  WGS  SOLUTION 

The  Unified  WGS  Solution  is  a  solution  for  geodetic  positions  and 
associated  parameters  based  on  an  optimum  combination  of  available  data. 
Initially,  a  normal  equation  matrix  was  formed  based  on  each  of  the 
previously  mentioned  data  sets.  Then,  the  individual  normal  equation 
matrices  were  combined  and  the  resultant  matrix  solved  to  obtain  the  WGS 
parameters.  In  order  to  merge  or  combine  the  individual  normal  equation 
matrices,  the  following  information  is  needed  or  required: 

-  The  normal  equations  themselves  free  of  extraneous 
cond i f i ons . 
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The  number  of  degrees  of  freedom. 


-  The  weighted  sum  of  the  squares  of  the  discrepancies 
(adjusted  if  necessary  for  any  eliminated  unknowns). 

-  An  estimate  of  the  variance  of  unit  weight. 

-  Initial  estimates  of  the  parameters. 

Given  this  information,  a  Unified  WGS  Solution  is  theoretically  very 
simple.  The  normal  equations  for  each  data  set  are  adjusted  so  that  the 
unknowns  are  homogeneous;  i.e.,  they  must  be  corrections  to  the  same 
initial  values  and  in  the  same  units.  Statistically,  they  are  homogenized 
by  dividing  each  by  its  variance  of  unit  weight.  Then,  it  is  only 
necessary  to  sum  the  various  normal  equation  matrices  and  solve  for  the 
unknowns . 

Unfortunately,  all  of  the  required  information  was  not  available  for 
all  the  normal  equation  matrices.  In  such  cases,  the  only  recourse,  with¬ 
out  going  to  the  expense  of  reforming  the  normal  equations,  was  to  adjust 
the  weights  assigned  to  each  system  to  bring  the  results  into  conformity 
with  reality.  That  is,  various  unified  solutions  were  made  using  differ¬ 
ent  weights  and  the  results  were  tested  against  external  standards.  As 
one  test,  for  example,  the  geopotential  coefficients  from  unified  solutions 
were  used  to  determine  goodness-of-fit  with  respect  to  ground-based 
satellite  tracking  data  for  an  actual  orbit.  Of  course,  each  unified 
solution  was  tested  for  consistency  with  the  individual  data  sets  which 
it  incorporated. 

Most  of  the  data  sets  are  complementary  to  each  other,  but  limited 
in  some  way  in  their  ability  to  uniquely  define  WGS  parameters.  For 
example,  due  to  the  lack  of  high  quality  worldwide  gravity  coverage,  in 
particular  in  the  southern  hemisphere,  existing  mean  free  air  gravity 
anomaly  data  poorly  defines  the  zonal  and  lower  degree  tesseral  harmonic 
coefficients  of  the  geopotential.  Similarly,  many  of  the  higher  degree 
tesseral  harmonic  coefficients  are  not  well  defined  from  satellite 
tracking  data  due  to  the  small  magnitude  of  the  orbital  perturbations 
at  geodetic  satellite  altitudes.  In  addition,  a  sufficient  number  of 
satellites  at  distinctly  different  inclinations  is  not  available  to  reduce 
the  correlation  between  the  coefficients.  Factors  such  as  these  provide 
the  logic  which  justifies  a  Unified  WGS  Solution  which  exploits  the 
advantages  and  at  the  same  time  minimizes  the  disadvantages  of  each  data 
set. 


Data  sets  available  in  the  form  of  normal  equations  matrices  and  in 
addition  to  ground  survey  information,  were  Doppler,  surface  gravity, 
astrogeodetic,  SECOR,  BC-4,  and  SAO  (Standard  Earth  II).  The  surface 
gravity  matrix  had  only  gravitational  model  parameters.  The  astro- 
geodetic  matrix,  in  addition  to  gravitational  model  parameters,  contained 
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unknowns  representing  datum  shift  parameters  and  a  change  in  the  semi¬ 
major  axis.  SECOR  and  BC-4  had  only  satellite  tracking  station  parameters, 
and  Doppler  and  Standard  Earth  II  included  both  satellite  tracking  station 
and  gravitational  model  parameters.  The  ground  survey  data  was  used  to 
provide  a  tie  between  stations  of  the  different  data  sets;  for  example, 
Doppler  and  BC-4,  Doppler  and  SECOR,  and  SECOR  and  BC-4.  However,  to 
minimize  the  effect  of  possible  survey  errors  on  the  solution,  the  use  of 
ground  survey  information  to  form  condition  equations  was  restricted  to 
collocated  stations.  This  does  not  eliminate  the  use  of  condition  equations 
for  precisely  surveyed  long  baselines  to  better  establish  scale  for  the 
system. 

A  starting  assumption  for  the  Unified  WGS  72  Solution  was  that  the 
normal  equations  obtained  from  local  survey  data  and  geometrically  pro¬ 
cessed  satellite  observational  material  (BC-4  and  SECOR),  could  be  taken 
at  face  value.  Due  to  neglect  of  gravitational  and  other  model  errors, 
normal  equations  obtained  from  dynamic  satellite  methods  and  geophysical 
data  reduction  might  require  some  deweighting.  For  the  surface  gravity 
data,  modification  of  data  weights  to  include  a  spatially  correlated 
model  term  with  appropriate  standard  error  and  decorrelation  distance 
would  have  given  more  balance  to  the  worldwide  data  representation  by 
that  matrix.  Because  the  astrogeodetic  normal  equations  were  not  formed 
from  the  measurements  themselves,  but  from  equally  weighted  observations 
of  the  finished  astrogeodetic  geoid  at  selected  locations,  important 
strengths,  weaknesses,  and  singularities  of  the  data  have  been  obscured. 

In  examining  the  matrices  used,  it  was  noted  that  only  the  SAO 
Standard  Earth  II  and  Doppler  matrices  required  joint  improvement  of  both 
stations  and  gravitational  parameters.  (The  connection  between  astro¬ 
geodetic  datum  shifts  and  station  positions  was  never  completed  with  a 
datum  survey  normal  matrix.)  Experience  with  Doppler  matrices  has  shown 
that  only  limited  correlation  exists  between  gravitational  coefficient 
and  station  coordinate  improvements.  Further,  the  Doppler  matrix  was 
expected  to  be  a  strong  contributor  to  the  unified  solution.  Therefore, 
the  basic  procedure  adopted  was  a  separate  determination  of  the  weights 
of  station-only  or  gravity-only  matrices  in  preliminary  solutions, 
followed  by  the  full  combined  solution  for  all  parameters.  In  the 
preliminary  phase,  the  SAO  Standard  Earth  II  matrix  was  treated  with  the 
other  gravity  determining  >■  "rices. 

As  weights  are  changed  in  combined  solutions,  each  set  of  normal 
equations  prefers  to  yield,  along  its  own  ill-determined  parameter, 
directions  which  can  make  the  unit  variances  insensitive  acceptability 
criteria.  As  alternative  criteria  for  judging  station  solutions,  the  RMS 
deviations  of  the  trial  solution  coordinates  were  expected  not  to  exceed 
similar  deviations  of  Doppler  coordinates  (adjusted  by  survey  differences 
to  the  non-Doppler  positions)  from  the  originators  sets.  According  to 
this  standard,  the  highest  satisfactory  trial  weight  for  the  BC-4  (and 
survey)  matrix  was  10,  relative  to  the  Doppler  weight.  The  SECOR  equations 
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were  entered  into  this  combination  at  the  highest  weight  relative  to  BC-4 
recommended  by  DMATC  and,  as  disturbances  were  minimal,  this  completed 
the  relative  weighting  determination  for  station  matrices.  Table  2 
compares  residuals  achieved  with  the  goals  for  each  data  set. 


Table  2 


Residuals  From  Station  Solutions 


Data  Set 

Type  of  Residual 

Goal 

Achieved 

Doppler 

2.2  m 

2.2  m 

BC-4 

RMS  Deviation  From 

6.6 

5.4 

SECOR 

Originators  Solution 

12.0 

10.0 

SAO 

9.0 

9.0 

Su  rvey 

■  _  —  ..  ... - 

RMS  Over  54  Sets 

Of  Tied  Stations 

0.5  m 

0.6  m 

Scalar 

RMS  of  8  Lines 

2  x  10"6 

1.5  x  10-6 

The  calibration  consisted  of  a  determination  of  biases  for  scale 
and  longitude  for  transformation  of  NWL-9D  positions  to  WGS  72.  Originally, 
the  scale  correction  was  based  on  the  comparison  of  satellite  and 
terrestrial  results  as  previously  discussed.  Since  the  equatorial 
radius  corresponding  to  the  dynamic  geoid  defined  by  an  interim  solution 
(WGSN-44)  geopotential  coefficients  and  the  mean  sea  level  values  of 
selected  Doppler  stations  is  6378133  meters,  a  scale  correction  is 
necessary  for  conversion  to  the  adopted  WGS  72  semimajor  axis  (6378135 
meters).  The  longitude  correction  was  defined  by  enforcing  a  selected 
set  of  Doppler  coordinates  to  agree  with  gravimetrical ly  derived  values 
in  the  NAD  27  area. 

The  merger  was  effected  through  an  adjustment  of  scale,  longitude, 
and  z-axis  origin  of  WGSN-44,  and  justified  by  the  low  residual  difference 
(2.2  meters,  RMS)  between  the  WGS  72  and  WGSN-44  (adjusted)  station 
coordinates  used  in  determining  the  transformation  constants.  Of  these, 
the  scale  and  z-shift  were  required  to  compensate  biases  inherited  from 
the  Doppler  normal  equation  matrix,  while  the  longitude  adjustment  was 
necessary  to  allow  for  the  difference  between  the  gravimetric  and  satellite 
standards.  Table  3  contains  the  formulas  and  parameters  for  relating 
NWL-9D  coordinates  to  WGS  72. 


A 


Table  3 

Formulas  and  Parameters 
To  Transform  NWL-9D  Coordinates 
To  WGS  72  Coordinates 


Formulas 

A<t>"  =  Af  sin  2$/s1n  1" 

AX"  =  0V260 

aH  -  a  Af  s  i  n2 (f>  -  Aa  +  Ar 

Pa rameters 

Ar  =  -5.27  meters 

Aa  =  -10.0  meters 

Af  =  -0.112415  x  10“6 

a  =  6378135  meters 

Instructions 

To  obtain  WGS  72  coordinates,  add  the 

A<t>,  AX,  aH  corrections  to  the  NWL-9D 
coordinates  (<t>,  x,  and  H,  respectively). 
Latitude  is  positive  North  and 

Longitude  is  positive  East. 

\ 
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MGS  72  PARAMETERS 


Mean  Earth  Ellipsoid 


In  determining  the  MGS  72  Ellipsoid  and  associated  parameters,  the 
Committee  decided  quite  early  to  closely  adhere  to  the  thoughts  and 
approach  used  by  the  International  Union  of  Geodesy  and  Geophysics  (IUGG) 
in  establishing  the  Geodetic  Reference  System  1967  (GRS  67)  [11].  Accord¬ 
ingly,  an  equipotential  ellipsoid  of  revolution  was  taken  as  the  form  for 
the  MGS  72  Ellipsoid.  An  equipotential  ellipsoid  is  simply  an  ellipsoid 
defined  to  be  an  equipotential  surface;  i.e.,  a  surface  on  which  all  values 
of  the  potential  are  equal.  Given  an  ellipsoid  of  revolution,  it  can  be 
made  an  equipotential  surface  of  a  certain  potential  function,  the  normal 
gravity  potential,  U.  This  normal  gravity  potential  can  be  uniquely 
determined,  independent  of  the  density  distribution  within  the  ellipsoid, 
by  using  any  system  of  four  independent  parameters  as  the  defining  con¬ 
stants  of  the  ellipsoid.  To  determine  the  normal  gravity  potential  with¬ 
out  resorting  to  the  ;se  of  a  mass  distribution  model  for  the  ellipsoid, 

U  can  be  expanded  into  a  series  of  zonal  ellipsoidal  harmonics  of  linear 
eccentricity  in  (a2-b2)5j.  The  coefficients  in  the  series  are  determined 
by  using  the  condition  that  the  ellipsoid  is  an  equipotential  surface 
(U  =  constant).  Since  all  the  zonal  coefficients  vanish,  except  the  two 
of  degree  zero  and  two,  a  closed  finite  expression  is  obtained  for  U. 

Normal  gravity  (y),  the  gradient  of  U,  is  given  at  the  surface  of  the 
ellipsoid  by  the  closed  formula  of  Somigliana  [12]: 


ave  cos2#  +  byp  sin2# 
(a2cos2#  +  b2  sin2#)1* 


where 

a  =  semimajor  axis  of  the  ellipsoid 
b  =  semiminor  axis  of  the  ellipsoid 
Yfi  *  normal  gravity  at  the  equator 

Yp  *  normal  gravity  at  the  poles 

#  =  geodetic  latitude 

Thus,  the  equipotential  ellipsoid  serves  not  only  as  the  reference  surface 
or  geometric  figure  of  the  earth  but  leads  to  a  closed  formula  for  normal 
gravity  at  the  ellipsoidal  surface,  a  formula  easily  modified  for  spatial 
applications. 
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Consistent  with  the  IUGG  definition  of  GRS  67,  the  Committee  took  the 
four  defining  parameters  of  the  WGS  72  Ellipsoid  to  be  the  semimajor  axis 
(a),  the  earth's  gravitational  constant  (GM)  and  angular  velocity  (u),  and 
the  second  degree  zonal  harmonic  coefficient  of  the  geopotential  ( C ^  Q). 

Other  parameters  associated  with  the  ellipsoid,  such  as  the  semiminor 
axis  (b)  and  the  flattening  (f),  including  the  normal  gravity  formula,  are 
calculated  using  the  defining  parameters.  These  and  other  parameters 
associated  with  the  ellipsoid  are  given  in  Table  4. 

Semimajor  Axis 


The  value  adopted  by  the  Committee  for  the  semimajor  axis  (a)  of  the 
WGS  72  Ellipsoid  is 

a  =  6378135  ±5  meters. 

The  adoption  of  an  a-value  10  meters  smaller  than  that  for  the  WGS  66 
Ellipsoid  was  based  on  several  calculations  and  indicators.  One  of  the 
more  extensive  calculations  was  made  using  an  approach  based  on  the 
general  combination  of  satellite  and  surface  gravity  data  for  position 
and  gravitational  field  determinations  [13].  There  are  two  variations  of 
the  general  combination  procedure  differing  principally  in  that  one  does 
not  involve  a  gravitational  field  determination.  Using  that  procedure  and  a 
related  computer  program,  various  sets  of  satellite  derived  station 
coordinates  and  gravimetric  deflection  of  the  vertical  and  geoid  height 
data  were  used  to  determine  local -to-geocentric  datum  shifts,  datum 
rotation  parameters,  a  datum  scale  parameter  and  a  value  for  the  semimajor 
axis  of  the  WGS  Ellipsoid.  Eight  solutions  were  made  with  the  various 
sets  of  input  data,  both  from  an  investigative  point  of  view  and  also 
because  of  the  limited  number  of  unknowns  which  could  be  solved  for  in  any 
individual  solution  due  to  computer  program  limitations.  Doppler  satellite 
tracking  and  astro-gravimetric  datum  orientation  stations  were  used  on  a 
selected  basis  in  the  various  solutions.  In  these  eight  solutions,  the 
input  values  for  the  semimajor  axis  and  flattening  of  the  ellipsoid  were 
6378145  ±15m  and  1/298.26,  respectively.  Different  combinations  of  Doppler 
satellite  tracking  and  astro-gravimetric  datum  orientation  stations  pro¬ 
vided  values  for  the  semimajor  axis  ranging  from  6378134.5  to  6378137.2 
meters.  Also,  eight  additional  solutions  were  made  in  which  the  only 
change  was  in  the  input  value  for  the  ellipsoidal  semimajor  axis.  Using 
6378130  ±15m  for  this  input  parameter,  the  solutions  for  the  semimajor  axis 
ranged  from  6378133.6  to  6378136.6  meters.  Based  on  these  results  and 
those  from  other  related  studies  accomplished  by  the  Committee,  the  value 
a  =  6378135  meters  was  adopted. 
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Table  4 


WGS  72  Ellipsoid 


Geodetic  and  Geophysical  Parameters 

Parameters 

| 

Notation 

Magis  i  tude 

Gravitational  Constant* 

GM 

398600.5  km3/sec2 

±0.4 

Second  Degree  Zonal* 

^2,0 

-484. 1605xl0”6 

-- 

Angular  Velocity* 

U> 

0.7292115147x10  4 
rad/ sec 

±0.1x10” 13 

Semimajor  Axis* 

a 

6378135  meters 

±5 

Flattening 

f 

1/298.26 

±0.6x10” 7 

Equatorial  Gravity 
(Absolute  System) 

Ye 

978033.26  mgal 

±1.8 

Gravitational  Constant 
(Mass  of  Earth's 
Atmosphere  Included) 

GM‘ 

398600.8  km3/sec2 

±0.4 

♦The  defining  parameters  of  the  WGS  72  Ellipsoid. 
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Table  4  (Cont'd) 
WGS  72  Ellipsoid 


Associated  Constants 

Constants 

Notation 

Formula 

Value 

Semiminor  Axis 

b 

■ 

1 

6,356,750.5  m 

Major  Eccentricity 

e 

e  =  [^(2-^)3*= 

0.08181881066 

Major  Eccentricity 
Squared 

e2 

e2  =  f (2-f ) 

0.006694317778 

Minor  Eccentricity 

e' 

■BBS 

0.08209405392 

— 

Axis  Ratio 

b/a 

b/a  =  1-f 

0.9966472205 

Radius  of  Sphere  with  Equal  Area  Rft 

6,371,005.2  m 

Radius  of  Sphere  wi 

in  Equa  ¥ 

n 

ol  ume  Rv 

6,370,998.9  m 

Ellipsoid  Potential 

Uo 

6,263,688  kgal  m 

Mean  Value  of  Norm*.  Gravity  ^Potsdam 

979,772.87  mgal 

Mean  Value  of  Normal  Gravity  ^Absolute 

979,758.87  mgal 
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Gravitational  Constant 


The  earth's  gravitational  constant  (GM),  the  product  of  the  universal 
gravitational  constant  and  the  earth's  mass,  can  be  obtained  by 
different  techniques  [14].  Ouring  recent  years,  the  most  reliable 
values  of  GM  have  resulted  from  the  analysis  of  lunar  and  planetary 
probe  tracking  data.  For  example,  the  GM  value  determined  by  the 
Jet  Propulsion  Laboratory  (JPL)  from  the  combination  and  analysis  of 
radio  tracking  data  from  four  lunar  missions.  Rangers  VI  -  IX,  was 
adopted  for  use  with  WGS  66.  Of  the  similar  more  recent  determinations, 
JPL  has  identified  the  GM  value  determined  from  Mariner  9  radio  tracking 
data  as  probably  a  more  accurate  result  to  date  [15].  This  value 

GM  =  398600.8  +0.4  km3/sec; 

which  contains  the  mass  of  the  earth's  atmosphere,  is  one  of  the  two  GM 
values  adopted  for  use  with  WGS  72.  The  second  WGS  72  GM  value  does_not 
include  the  mass  of  the  earth’s  atmosphere  and,  along  with  a,  u  and  q, 

is  one  of  the  four  defining  parameters  of  the  WGS  72  Ellipsoid.  It  is 
discussed  later. 

In  determining  GM,  JPL  used  the  value 
c  =  299792.5  '0.3  km/sec 

for  the  velocity  of  light.  Since  that  time,  new  determinations  of  the 
value  for  c  have  been  reported  by  the  National  Bureau  of  Standards  (NBS) 
[16]  [17].  The  two  values  are: 

c  =  299792.4562  *0.0011  km/ sec 

c  =  299792.462  ±0.018  km/sec 

In  June  1973,  the  value  of  c  was  revised  by  the  Consultative  Committee 
for  the  Definition  of  the  Meter  (under  the  International  Bureau  of 
Weights  and  Measures).  This  action  occurred  as  a  result  of  their 
selection  of  a  recommended  value  for  the  wavelength  of  the  methane- 
stabilized  laser,  which  thereby  (in  conjunction  with  NBS  laser  frequency 
measurements)  established  a  recommended  value  for  c  [18].  The  recommended 
value  was  adopted  by  the  Committee  for  use  in  D0D  applications: 

c  =  299792.458  *0.0012  km/sec 

As  previously  stated,  the  above  JPL  determined  value  for  GM  contains  the 
mass  of  the  earth's  atmosphere.  To  derive  a  GM  value  which  does  not 
contain  the  mass  of  the  earth's  atmosphere,  values  are  needed  for  the 
universal  gravitational  constant  (G)  and  the  mass  of  the  earth's 
atmosphere  (ma)-  The  value 
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m  =  (5.136  +0.007)  x  1021  grams 

d 

was  selected  along  with 


G  =  66.720  +0.041  x  10'9  cm3/sec2 

for  this  purpose.  The  rn  value,  above,  was  selected  on  the  basis  of  the 

a 

greater  rigor  applied  in  its  determination.  The  product  of  G  and  m  ,  when 

d 

subtracted  from  GM,  gives 

GM  =  398600.5  +0.4  km3/sec2 

The  "esult  is  a  GM  value  which  does  not  contain  the  mass  of  the  earth's 
atmosphere  and  is  also  adopted  for  use  with  WGS  72. 

The  selection  of  the  proper  value  depends  on  the  type  of  application 
considered.  For  example,  some  satellite  and  space  studies  require  the 
GM  which  includes  the  earth's  atmosphere.  In  geodetic  computations 
involving  the  normal  potential  (gravity  formula,  gravity  anomalies, 
geoid  heights,  etc.)  it  is  proper  to  use  the  GM  excluding  atmosphere. 

Angular  Velocity 

The  angular  velocity  of  the  WGS  72  Ellipsoid  was  taken  to  be  equal  to 
that  of  the  actual  earth.  The  value  selected  for  the  WGS  72  Ellipsoid 

w  =  7.292115147  x  10~5  radians/second 

is  a  rounded  value  of  the  GRS  67  value  which  has  an  extra  digit  [11]. 

The  GRS  67  value  was  obtained  by  two  solutions  and  is  consistent 
with  the  International  Astronomical  Union  (IAU)  System  of  Astronomical 
Constants,  which  were  adopted  by  the  IAU  in  1964  [19]. 

Second  Degree  Zonal  Harmonic  Coefficients  of  the  Geopotential 

_  The  second  degree  zonal  harmonic  coefficient  of  the  geopotential  is  the 
C^  g  value  in  the  following  series  where  the  geopotential  form  is 

restricted  to  zonal  harmonics: 


U  = 


GM 


1  + 


l 

n=2 


Cn,0  (f) 


n  _ 


Vo  (s,n*') 


(2) 


Due  to  the  limitation  of  available  surface  gravity  data,  satellite 
observational  data  provides  the  best  means  for  obtaining  coefficients 
for  equation  (2).  The  Cn  g  coefficients  are  expressed  in  terms  of  the 


rates  of  change  of  the  orbital  elements,  the  changes  in  the  elements 
being  observed  after  several  revolutions  rather  than  directly.  Values 
for  the  even  degree  coefficients  are  principally  influenced  by  secular 
or  long-periodic  variations  in  the  ascending  node  and  perigee,  odd  degree 
coefficients  are  principally  influenced  by  similar-type  variations  in  the 
eccentricity  and  inclination. 

Using  satellite  data  and  an  equation  such  as  equation  (2),  various 
determinations  have  been  made  of  a  finite  number  of  zonal  harmonic 
coefficients.  For  example,  King-Hele  [20]  [21]  and  Kozai  [22]  made 
several  determinations  during  the  mid  and  late  1960s.  Also,  using 
the  more  general  expression  for  the  geopotential. 


J2  JQ  <?>"  P™  (C™  cos  m>  *  snm  s1n  mi> 


(3) 


the  zonal  harmonic  coefficients  are  determined  along  with  the  tesserals  in 
the  standard  least  squares^  solution  for  the  gravitational  model  coefficients 
(^nm*  snm)*  Selectin9  a  C2  0  value  to  serve  as  one  the  four  defining 

parameters  of  the  WGS  72  Ellipsoid  is  more  desirable  from  a  consistency 
standpoint  since  it  is  part  of  a  larger  set  of  correlated  coefficients 
forming  the  earth's  gravitational  model.  A  C2  g  value  selected  from 

a  solution  involving  equation  (2)  would  not  necessarily  be  consistent 
with  the  coefficient  set  from  a  solution  involving  equation  (3). 

Therefore,  theCommi ttee  adopted  as  a  defining  parameter  of  the  WGS  72 
Ellipsoid  the  C2  Q  value  from  the  geopotential  coefficient  set  selected 

as  the  WGS  72  Earth  Gravitational  Model: 

C2>0  =  -484.1605  x  10'6 

Flattening  (Ell ipticity) 

The  flattening  (f)  or  ell ipticity  of  the  ellipsoid  is  related  to  the 
second  degree  zonal  harmonic  coefficient  of  the  geopotential  by  the 
equation 


C  =  -  midJ/im  +  !PlH3/Um-t2/7lfl  +  Q(f3) 

2,0  3  3 

where  m  =  to2  a/ye.  Using  WGS  72  values  for  C2  Q  and  each  parameter  in 

the  equation  for  m,  equation  (4)  was  solved  for  f  by  iteration.  From 
this  computation,  the  flattening  of  the  WGS  72  Ellipsoid  is: 
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f  =  1/298.2638 


This  value  has  been  rounded  to  1/298.26  which  is  the  adopted  value. 

Normal  (Ellipsoid al )  Gravity  Formula 

The  WGS  72  Ellipsoidal  Gravity  Formula  is  in  the  form  of  the  Chebychev 
approximation  [11]: 

Y  =  A  (1  +  c  sin2$  +  c4  sin4<j>)  (5) 

where 

A  =  ye  ^  +  32  a<  ^ 


c2  ~  a?  "  ae 


c4  =  a4  +  j  a6 

The  coefficients  for  this  formula  were  obtained  by  using  the  WGS  72 
parameters,  resulting  in  the  WGS  72  Gravity  Formula: 

Y  *  978.03327  (1  +  0.005278994  sin2<j>  +  0.000023461  sin4$)  gal 
where  1  gal  =  1  cm/sec2  and  1  milligal  (mgal)  =  10  gal. 

The  WGS  72  Geo  id 

In  geodetic  appl ications,  three  different  surfaces  or  earth  figures 
are  normally  invoked.  In  addition  to  the  earth's  natural  or  physical 
surface,  these  include  a  geometric  or  mathematical  reference  surface, 
the  ellipsoid,  and  an  equipotential  surface  called  the  geoid.  The 
geoid  theoretically  coincides  with  mean  sea  level  over  the  oceans  and 
extends  hypothetically  beneath  all  land  surfaces.  Normally,  the  geoid 
makes  its  appearance  in  geodetic  positioning  applications  by  the 
relationship 

H  =  N  +  h  (6) 


where 


H  =  geodetic  height 

N  =  geoid  height 

h  =  height  ubnve  mean  sea  level. 
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For  general  purposes,  the  form  of  the  geoid  is  usually  indicated  by 
a  contour  chart  which  shows  the  location  of  the  geoid  above  or  below  a 
geometric  reference  surface  or  figure  (the  ellipsoid)  selected  to  approxi¬ 
mate  the  earth.  When  available,  an  astrogeodetic  contour  chart  is  the 
first  practical  choice  used  to  obtain  geoid  height  data  (Figure  5). 
Although  astrogeodetic  charts  are  very  detailed  and  accurate,  there  are 
gaps  in  the  coverage.  For  practical  applications,  geoid  heights  can  be 
computed  mathematically  from  the  appropriate  formulation  rather  than 
scaled  from  a  worldwide  geoid  contour  chart. 

THE  CONVERSION  OF  GEODETIC  DATUM  TO  WGS  72 

In  the  development  of  local-to  WGS  72  datum  shifts,  different  geo¬ 
detic  disciplines  have  been  used  and  the  results  from  each  investigated, 
analyzed  and  compared.  This  redundancy  of  disciplines  and  data  provides 
assurance  that  the  system  accepted  as  WGS  72  is  the  best  attainable  using 
the  methods  and  data  available  in  1972.  Noting  the  variations  in  results 
obtained  from  the  different  disciplines  was  helpful  in  assigning  accuracy 
statements  to  the  adopted  parameters.  Table  5  compares  results  from 
gravimetric  data,  Doppler  satellite  tracking  data,  and  a  combination  of 
the  two  sources  used  to  orient  four  principal  datums  to  a  geocentric 
system.  (The  local-to-WGS  66  datum  shifts  are  shown  for  comparison 
purposes.)  Corrections  to  shift  the  previously  existing  operattonal 
Doppler  coordinate  set  (NWL-9D)  to  WGS  72  were  shown  in  Table  3.  Using 
these  corrections,  a  large  number  of  Doppler  TRANET  and  Geoceiver 
station  coordinates  was  available  worldwide  for  deriving  local-to-WGS 
72  datum  shifts.  The  derivation  of  the  shifts  is  discussed  below. 

NAD  27  to  WGS  72 


Stations  distributed  over  the  continental  United  States  and  Mexico 
were  selected  from  the  available  coordinates.  The  selection,  listed  in 
Table  6,  consists  principally  of  Geoceiver  stations  with  some  permanent 
Doppler  Tracking  Network  (TRANET)  stations  and  ITT-5500  stations.  Four 
locations  (Howard  County,  MD;  Beltsville,  MD;  Las  Cruces,  NM;  and  Moses 
Lake,  WA)  have  both  TRANET  and  Geoceiver  stations.  These  stations  were 
used  to  derive  datum  shifts  for  the  NAD  27  area,  excluding  Alaska  and 
Canada.  Figures  6,  7  and  8  show  the  NAD  27-to-WGS  72  datum  shifts  in 
terms  of  contours  of  the  aX,  aY,  aZ  components.  The  datum  shifts 
contained  in  Tables  5  and  7  are  mean  values  of  the  datum  shifts  for 
individual  stations.  Variations  of  several  meters  (up  to  10  m)  exist 
between  the  mean  NAD  27-to-WGS  72  datum  shifts  and  the  shifts  observed 
at  each  location.  The  standard  deviations  from  the  mean  are  o ^  =  ±4.4  m, 

oaY  =  ±3.5  m,  o  ^  -  *3.9  m.  Thus,  while  the  mean  NAD  27-to-WGS  72  datum 
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Figure  5.  Astrogeodetic  WGS  72  Geoid  (meters) 


^Longitude  constraint  imposed 


Table  6 


Stations  Used  to  Derive  WGS  72  Datum  Shifts 


NAD  27  Area,  Alaska  and  Canada  Excluded 


Station 

Number 

Location 

Station 

Number 

Location 

103 

ill 

311 

321 

738 

742 

10003 

10006 

10018 

10019 

10020 
10021 
10023 
10029 
10031 
10046 
10055 
20000 
20001 

Las  Cruces,  NM 

Howard  Co.,  MD 
Prospect  Harbor,  ME 
Rosemount,  MN 

Moses  Lake,  WA 
Beltsville,  MD 
Greenville,  MS 

Tipton,  KS 

Jonestown,  TX 

Frankton,  IN 
Marysville,  IN 

Summit,  KY 

Mathiston,  MS 

Patrick  AFB,  FL 
Goldstone,  CA 

Salt,  NM 

Pillar  Point,  CA 
Howard  Co.,  MD 
Beltsville,  MO 

20002 

20003 

20015 

20016 
30025 

30027 

30028 

30029 

30030 

30031 

30032 

30033 

30034 
30038 
30078 
30089 

30098 

30099 
8004 

_ 

Las  Cruces,  NM 
Wrightwood,  CA 
Woodbine,  GA 

Columbia,  MS 
Bloomfield,  OH 
Greenville,  OH 
Metamora,  IL 

Moses  Lake,  WA 

Green  River,  UT 
Seligman,  AZ 

Navajo,  AZ 

Miami,  AZ 

Ajo,  AZ 

Blaine,  WA 

Westford,  MA 

Los  Mochis,  Mexico 
Orland,  CA 

Malta,  MT 

Marshall  Point,  ME 

NAD  27  Area,  Alaska  and  Canada  Only 

Station 

Number 

Location 

Station 

Number 

Location 

14 

740 

807 

30050 

30051 

30052 

Anchorage,  AK 

Nome,  AK 

Hall  Beach,  Canada 
Kotzebue,  AK 

Clear,  AK 

Point  Barrow,  AK 

30053 

30054 
30056 
30065 
30067 

_ 

Barter  I . ,  AK 

Boundary,  AK 

Cold  Bay,  AK 

Cape  Sarichef,  AK 
Gilmore  Creek,  AK 

26 
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Figure  6.  Contour  Chart  -  NAD  27  to  WGS  72  Datum  Shifts,  Component  (Contour  Interval  =  2  meters). 
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Contour  Chart  -  NAD  27  to  WGS  72  Datum  Shifts,  Component  (Contour  Interval  =  2  meters). 


Table  7 

Datum  Shift  Constants 
{Geodetic  Datum  to  WGS  72) 


Geodetic  Datums 
and 

Constants 

_ 

Reference  Ellipsoids 

wm\ 

mm 

EJBI 

AfXlO"4 

North  American  1927 

-  22* 

157* 

176* 

-  71.400 

-0.37295850 

{Clarke  1866) 

Alaska  and  Canada 

-  9 

139 

173 

-  71.400 

-0.37295850 

European 

(International ) 

-  84 

-103 

-127 

-253.000 

-0.14223913 

Tokyo 

(Bessel ) 

-140 

516 

673 

737.845 

0.10006272 

Australian  Geodetic 
(Australian  National) 

-122 

-  41 

146 

-  25.000 

-0.00112415 

Ordnance  Survey  of  Great 
Britain  1936  (Airy) 

368 

-120 

425 

571.604 

0.11928812 

South  American  1969 
(South  American  1969) 

-  77 

3 

-  45 

-  25.000 

-0.00112415 

Old  Hawaiian 

(Clarke  1866) 

Maui 

65 

-272 

-197 

-  71.400 

-0.37295850 

Oahu 

56 

-268 

-187 

-  71.400 

-0.37295850 

Kauai 

46 

-271 

-181 

-  71.400 

-0.37295850 

Johnston  Island  Astro  1961 
(International ) 

192 

-  59 

-211 

-253.000 

-0.14223913 

Wake-Eniwetok  1960  (Hough) 

Kwajalein  Atoll 

68 

-  44 

-135.000 

-0.14223913 

Wake  Island 

62 

-  22 

-135.000 

-0.14223913 

Eniwetok  Atoll 

62 

-  38 

-135.000 

-0.14223913 

♦Mean  value  for  the  NAD  27  area  excluding  Alaska  and  Canada;  see  also 
Figures  6,  7,  and  8. 
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Table  7  (Cont'd) 

Datum  Shift  Constants 
(Geodetic  Datum  to  WGS  72) 


Geodetic  Datums 
and 

Constants 

Reference  Ellipsoids 

daoi 

WBBM 

Wake  Island  Astro  1952 
(International ) 

283 

-  44 

141 

-253.000 

-0.14223913 

Canton  Island  Astro  1966 
( International ) 

294 

-288 

-382 

-253.000 

-0.14223913 

Guam  1963 

(Clarke  1866) 

-  89 

-235 

254 

-  71.400 

-0.37295850 

Ascension  Island  Astro  1958 
( International ) 

-214 

91 

48 

-253.000 

-0.14223913 

South  Asia 

(Fischer  1960) 

21 

-  61 

-  15 

-  20.000 

0.00449585 

Nanking  1960 

(International) 

-131 

-347 

B 

-253.000 

-0.14223913 

Arc  1950 

(Clarke  1880) 

-129 

-131 

-282 

-114.145 

-0.54781925 

Adindan 

(Clarke  1880) 

-152 

to 

CM 

1 

212 

-114.145 

-0.54781925 

Mercury  1960 
(Fischer  1960) 

NAD  27  Area 

-  25 

46 

-  49 

m 

0.00449585 

ED  Area 

-  13 

-  88 

-  5 

0.00449585 

TD  Area 

18 

-132 

60 

-  31.0 

0.00449585 

Modified  Mercury  1968 
(Fischer  1968) 

NAD  27  Area 

-  4 

12 

-  7 

-  15.0 

0.00449585 

ED  Area 

-  3 

1 

-  6 

-  15.0 

0.00449585 

TD  Area 

22 

34 

2 

-  15.0 
_ 

0.00449585 
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shifts  are  valid  within  ±5  m  in  each  component,  users  may  prefer  the 
greater  precision  attainable  by  using  the  datum  shifts  shown  in  Figures 
6,  7,  and  8.  Such  shifts  are  more  closely  related  to  the  geodetic 
control  in  a  localized  area  in  that  survey  errors  which  accumulate  over 
long  distances  from  the  datum  origin  are  eliminated.  The  result  is  a 
form  of  "local"  to  WGS  72  datum  shift.  In  this  sense,  individual  NAD 
27  to  WGS  72  datum  shifts  are  also  shown  in  Table  7  for  Alaska  and 
Canada. 

ED  to  WGS  72 


The  number  of  Doppler  stations  for  ED  is,  of  course,  much  less 
than  that  available  in  the  NAD  27  area  and  the  distribution  less 
desirable.  A  total  of  eight  stations,  four  TRANET  and  four  Geoceiver 
stations,  was  selected  for  deriving  ED  to  WGS  72  datum  shifts.  Tables 
5  „>rd  7  contain  the  mean  of  the  aX,  aY,  aZ  shifts  for  the  eight  stations. 
The  mean  values  define  the  ED  to  WGS  72  datum  shifts. 

TD  to  WGS  72 


Datum  shifts  for  Tokyo  Datum  were  determined  in  a  manner  similar 
to  that  for  the  ED  area;  however,  only  four  stations  were  used  (two 
TRANET  and  two  Geoceiver).  The  TD  to  WGS  72  datum  shifts  are  shown 
in  Tables  5  and  7. 

AGD  to  WGS  72 


Four  TRANET  stations  were  used  to  derive  WGS  72  datum  shifts  for  the 
Australian  Geodetic  Datum  (AGD).  The  datum  shifts  are  contained  in  Tables 
5  and  7. 

Other  Datum  Conversions 


In  addition  to  the  datum  shift  constants  for  the  preceding  datums. 
Table  7  also  contains  constants  for  converting  other  significant  geo¬ 
detic  datums  to  WGS  72.  The  shifts  are  also  based  on  Doppler  derived 
coordinates  in  combination  with  conventional  survey  data.  Included  are 
datum  shifts  for  several  island  areas.  These  are  good  examples  of  the 
utilization  of  Geoceiver  surveys  to  accomplish  geocentric  positioning 
in  terms  of  the  WGS  72.  Several  smaller  datums  exist  for  which  there 
are  no  direct  ties  to  WGS  72.  Transformations  between  these  local 
datums  and  major  datums  have  been  derived  where  possible  using 
terrestrial  and  satellite  information.  Datum  shifts  for  this  purpose 
are  shown  in  Table  8.  For  information  purposes,  some  interdatum 
relationships  for  datums  in  Table  7  are  also  included. 
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Table  8 


Interdatum  Shifts  for  WGS  72 


Geodetic  Datums 

Constants 

From 

To 

AX(m) 

AY(m) 

AZ(m) 

Aa(m) 

AfXlO-4 

SAD  69 

NAD  27 

-  55 

-154 

-221 

46.4 

0.3718344 

Chua  Astro 

SAD  69 

-  77 

239 

5 

-228.0 

-0.1411150 

Yacare 

SAD  69 

-  90 

160 

78 

-228.0 

-0.1411150 

Campo 

Inchauspe 

SAD  69 

-  83 

130 

120 

-228.0 

-0.1411150 

Hito  XVIII 

SAD  69 

87 

198 

125 

-228.0 

-0.1411150 

Bogota 

SAD  69 

354 

288 

-283 

-228.0 

-0.1411150 

La  Canoa 

SAO  69 

-225 

102 

-326 

-228.0 

-0.1411150 

Aerodist 

SAD  69 

-222 

108 

-317 

-228.0 

-0.1411150 

Adindan 

ED 

-  68 

77 

339 

138.855 

-0.4055801 

Arc  1950 

Adindan 

23 

-105 

-494 

0 

0 

Arc  1950 

ED 

-  45 

-  28 

-155 

138.855 

-0.4055801 

South  Asia 

ED 

105 

42 

112 

233.0 

0.1467350 

Malayan  Revised 
Triangulation 

South  Asia 

-  33 

918 

30 

850.937 

0.2788057 
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Datum  Conversion  Formulas 


The  Molodensky  Datum  Conversion  Formulas  [23]  previously  selected 
for  use  with  WGS  66  from  an  analysis  of  several  datum  conversion  formulas 
are  retained  for  WGS  72  purposes.  To  accomplish  the  conversion,  the  data 
in  Table  7  and  the  reference  ellipsoid  constants  (Table  9)  are  used  in 
the  Molodensky  Datum  Transformation  Formulas  (Table  10).  In  the  selected 
rectangular  coordinate  system,  the  axes  are  defined  as  follows: 

X-axis  =  Intersection  of  the  plane  through  the  Greenwich  meridian 
and  the  plane  of  the  equator. 

Y-axis  =  Measured  90°  east  of  Greenwich  in  the  plane  of  the  equator. 

Z-axis  =  Coincident  with  the  earth's  axis  of  rotation. 

EVALUATION 

Accuracy  estimates  of  the  parameters  defining  WGS  72  are  included 
in  Table  4.  In  addition.  Figure  9  contains  the  uncertainties  of  the 
datum  shift  constants  for  selected  datums.  These  range  from  ±5  meters, 
one  sigma,  in  the  NAD  27  area  up  to  +.15  meters  in  the  TD  and  AGO  areas.  The 
accuracy  values  were  determined  in  the  WGS  72  development  primarily  by 
error  analyses  involving  Doppler  and  physical  geodesy  results.  A  relative 
evaluation  of  WGS  72  is  provided  by  a  comparison  with  other  contemporary 
geocentric  systems.  Table  11  compares  parameters  and  datum  shifts  between 
WGS  72  and  two  recent  solutions  each  of  the  Smithsonian  Astrophysical 
Observatory  (Standard  Earth  II,  1969  [7],  Standard  Earth  III,  1973  [24])  and 
the  Goddard  Space  Flight  Center  (GSFC  GEM  4,  1972  [25],  GSFC  1973  [26]). 
Generally,  the  datum  shifts  agree  well  in  the  NAD  27  area  with  larger 
differences  in  the  ED  and  TD  areas.  Additional  datums  are  compared  in 
Figure  10  between  WGS  72  and  GSFC  1973.  The  good  agreement  in  the  NAD 
27  area  includes  the  separate  shifts  for  Alaska  and  Canada.  With  some 
exceptions,  the  differences  are  generally  less  than  the  one  sigma 
accuracy  values  quoted  for  WGS  72  and  well  within  the  three  sigma  level. 

The  WGS  72,  adopted  by  the  DOD,  represents  a  state  of  the  art 
WGS  solution  using  worldwide  geodetic/gravimetric  surveys  and 
electronic/optical  satellite  data,  and  as  such,  should  be  considered 
a  candidate  system  for  use  in  the  readjustment  of  the  North  American 
triangulation.  Further,  the  WGS  72  associated  with  increased  worldwide 
use  of  the  geoceiver,  provides  a  unique  combination  which,  with  inter¬ 
national  cooperation,  could  rapidly  lead  to  a  true  world  datum. 
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Table  9 


Fischer  1968  6378150  1/298.3  0.335  232  986  9  1.003  363  606  ]  0.996  647  670  1  i  0.006  693  421  6  0.993  306  578  4 


Table  10 

COORDINATE  TRANSFORMATION  FORM  (’LAS 
OEOOETIC  OATI'M  TO  W<;s  72 


A.  The  Standard  Molodensky  Formulas 

A0"  -  {-AX  sin  <o  r-os  A  aY  sin  0  sin  A  r  aZ  cos  0 
4  Aa  (KNe2  sin  0  cos  </>)'a 

+  Af|RM(a/b)  +  Ks  1  b  a)j  sin  0  cos  0 j  •  |(KM  4  H)  sin  1']  1 

AA"  -  |  AX  sin  A  t  aY  cos  A|  .  |(RS  4  H)  cos  0  sin  l”)  1 

AH  =  aX  cos  0  cos  a  *  aY  cos  0  sin  A  4  aZ  sin  0 

-Aa  (a'RNi  4  Af  (h'a)  RN  sin2  0 

B.  The  Abridged  Molodcnsky  Formulas 

A0"  -  |  -AX  sin  0  cos  A  AY'  sin  0  sin  A  +  AZ  cos  0  4  (aAf  +.  fAal  sin  2  0] 

.  |R„  sin  l"|  1 

AA”  -  |  aX  sin  A  4  A\'  Cos  A |  .  (R^  cos  0  sin  1"|  1 

AH  -  A\  cos  0  cos  A  4  aY'  cos  0  sin  A  4  aZ  sin  0  *  (aAf  4  fAal  sin2  0  -  Aa 

('.  Definition  of  Terms  in  ihe  Molodcnsky  Formulas 

0,  A.H  •  geodetic  coordinates  (old  ellipsoid) 

0  -  geodetic  latitude.  The  angle  between  the  carih's  equatorial  plane  and  the  ellij>- 
soidal  normal  at  a  point  (measured  positive  meth  from  the  equator,  negative 
sot:'  !i). 

A  good-tie  longitude.  The  angle  la-tween  the  plane  of  the  'Greenwich  meridian  and 
(he  plane  of  the  geodetic  meridian  of  the  [sunt  t measure J  in  the  plane  of  the 
equator,  positive  oast  from  tlreenwichl. 

H  the  distance  of  a  [ioint  from  the  ellipsoid  measured  along  the  ellipsoidal  norma) 
through  the  point. 

I!  -  N  4*11 

‘Indicates  parameters  which  do  not  appear  in  the  Abridged  \|o|'nh  nsl<>  Formulas. 
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Table  10  (Cont'd) 


N  =  gemd-ellipsoid  separation.  The  distance  of  the  geoid  above 
I  h  Ni  or  below  i  N)  |  he  ellipsoid 

h  -  distance  ol  a  (Kuril  Ironi  Ihe  geoid  (elevation  above  or  lielmv 
mean  ea  ■>  h 


Ac-.  AA  Jill 


corrections  to  Uansform  the  geodetic  coordinates  from  the 
old  datum  t"  \\  t  IS 


A  \  A^  A/ 

a 

*1. 


shifts  ik  tween  i  lipsoid  renters  of  the  old  datum  and  W(iS. 
sonumajor  axis  of  the  old  ellipsoid, 
semiminor  axis  of  i he  old  ellipsoii:. 


♦b  a  1  f 

f  *  flattening  of  the  old  elli(isoid. 

Aa.  Af  differenres  between  the  parameters  of  the  old  ellipsoid  and 
the  \\  t  iS  ellipsoid  t  \\  t  iS  in  mils  old  i 

e  -  montririH 


e  ■-  Sf  f 

Ks  radius  of  curvature  m  the  prime  vertical 


“KN  a  ( I  e  sin  ,.--i 

Ky  radius  of  eurv  at  on- m  t  tie  meridian. 
**RV  al  1  e’  l  ( I  e  mu  ,.-  i 1  ’ 


\i  'I  K:  All  A-quant  dies  tire  formed  bv  subtracting  old  ellipsoid  values  from  VVtiS  ellipsoid  values. 


‘Indicates  parameters  which  do  not  appear  m  Ihe  Abridged  Molodenskv  Formulas. 

**For  desk  calculator  computations  involiiiip  commonly  used  ellipsoids,  these  values  are  given  in 
(attitude  Funetion  Tables.  i . ■  •  .  I.atpofle  Kurtrlion.  ('lark*1  ) Slid  Spheroid.  AMS 'I’M  No.  t>S,  1S."7. 
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Figure  9.  WQS  72  Datum  Shifts  and  Accuracies  for  Selected  Datums. 
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Figure  10.  Comparison  of  Datum  Shifts,  WGS  72  Minus  GSFC  1973  (meters). 
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